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Abstract We consider a self-similar fragmentation process in which the generic particle of size x 
is replaced at probability rate by its offspring made of smaller particles, where a is some positive 
parameter. The total of offspring sizes may be both larger or smaller than x with positive probability. 
We show that under certain conditions the typical size in the ensemble is of the order t^-^/°' and that 
the empirical distribution of sizes converges to a random limit which we characterise in terms of the 
reproduction law. 

1 Introduction 

We study the following continuous-time model of particle fragmentation. Each particle in ensemble is 
characterised by a positive quantity which we call size. The generic particle of size x lives a random 
exponentially distributed time with parameter cc", where a is some fixed real number. During the life- 
time the size does not vary and at the end the particle splits into fragments of sizes x£,j, where {^j} is 
independent of the lifetime of the particle and follows a given probability distribution called reproduction 
law. Each particle is autonomous, meaning that the splitting probability rate and the descendant fragment 
sizes depend only on the size of the particle and not on the history of this or other coexisting particles. 
We are interested in the case a > 0, when particles with smaller size tend to live longer. We refer to the 
proceedings |14| and the survey [jP for a number of examples arising in physics and chemistry. 

The idea of the model was suggested by Kolmogorov in [22 , and the first results are due to Filippov 
[T^ . Brennan and Durrett ^| rediscovered an instance of the model in the context of binary interval 
splitting. In a recent series of papers Bertoin [^lEIIHl introduced more involved fragmentation processes 
in which a particle may produce infinitely many generations within an arbitrary time period or may 
undergo a continuous size erosion; see also [3 El El 1211 EDI ESI for related examples. 

The research so far was mainly focussed on the conservative case ~ 1 when the total size is 

preserved by each splitting. It has been shown that the particles demonstrate quite a regular long-run 
behaviour: the typical size in the ensemble is of the order and a scaled empirical distribution of 

sizes converges to a nonrandom limit. Filippov |19| proved the convergence of empirical distributions in 
probability (see also [5]), while Brennan and Durrett jl8| showed convergence with probability one in the 
binary case. Baryshnikov and Gnedin El studied a sequential interval packing problem which may be 
seen as a binary instance of dissipative fragmentation with X^Cj ^ 1 and F(X]Cj < 1) > Oj and proved 
convergence of the mean measures associated with the empirical distributions. We mention that some 
special cases of these mathematical results have also appeared in the literature in physics, see e.g. 01123 
and references therein. 

Conservative or dissipative fragmentations can be treated both as continuous-time interval splitting 
schemes, similar to discrete-time random recursive constructions (as in [2H|), or as state-discretised pro- 
cesses with values in Kingman's partition structures d El El ^1 ■ These approaches fail completely 
if the reproduction law allows the possibility of size creation, when the total size of the offspring may 
exceed the size of the parent particle. Such 'improper fragmentations' are both physically plausible and 
useful in the situations where the generalised size models some nonadditive quantity like, e.g. surface 
energy by aerosols. 
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By allocating particle of size a: at — log x the fragmentation process can be seen as a branching random 
walk with location-dependent sojourn times. From this viewpoint, a constraint on the sum of sizes seems 
rather odd, which suggests that such a condition is not essential for the asymptotics. In this paper we 
argue that this intuition is indeed correct, in the sense that the convergence of properly scaled empirical 
distributions of sizes holds under fairly general assumptions on the reproduction law. Though we do 
require that the individual offspring sizes cannot exceed the parent size, there is no constraint on the 
total offspring size. A new feature appearing in the general nonconservative case is that the limit of scaled 
empirical distributions is not completely deterministic, rather involves a random factor which admits a 
characterisation by a distributional fixed-point equation. This phenomenon reminds us, of course, of 
strong limit theorems for branching random walks; we shall discuss the connection and the differences in 
Section O 

The rest of this work is organised as follows. Notation and basic assumptions are given in Section 
12.11 and Section 12.21 introduces the genealogical structure and the so-called intrinsic martingale which 
plays a fundamental role in branching processes. Then we compute the first moment of power sums and 
then determine its asymptotics using a contour integral. This yields the convergence of mean measures 
in Section [3. 31 An alternative approach based on a limit theorem of Brennan and Durrett is presented in 
Section 13.41 The main result of convergence of scaled empirical measures is proved in Section l4.ll Then 
we provide some examples in Section|Sl and finally, in Section|H| we sketch the extension of the preceding 
results to self-similar fragmentations with possibly infinite reproduction measure. 



2 Preliminaries 

2.1 Definitions and assumptions on the reproduction law 

The collection {^j}jeN of offspring sizes of a unit particle is identified with a random sequence of non- 
negative real numbers ranked in the decreasing order (by convention, — when j is larger than the 
number of children). We shall also view {CjljsN as a random set, defined formally as a counting random 
measure ^ 5(^. on ]0, 1] . Basically we require that 

ij e [0, 1] , E#{j : > 0} > 1 , P(ei = 1) < 1 . (1) 

Many features of the fragmentation process can be expressed in terms of the structural measure 

a(i?)=E#{j:Oei?}, SC]0,1], 

and its Mellin transform 

m = / x^a{Ax) = E ^ 

which we call the characteristic function. In particular, the conditions in JQ) amount to the assumptions 
that a is supported by [0, 1], that crjO, 1] > 1 and that ^{l} < 1. 

Because < (/)(5R/3), the natural domain of definition of </) is a complex halfplane to the right of 

the convergence abscissa Pa of the integral. If (3a = —00 the halfplane is the whole plane, and otherwise 
the halfplane may be open or closed. The characteristic function is analytical in the halfplane, strictly 
decreasing on the real axis, and in view of Q satisfies 0(0) > 1 and ^ o-{l} < 1 as 00. 

It is crucial for our results and will be assumed throughout that there exists the Malthusian exponent 
(3* > satisfying the equation 

m - 1 . (2) 

If the Malthusian exponent exists then it is unique and there are no solutions to ^ in the halfplane 
> (3*. And if some (3 ^ (3* with = (3* satisfies (jJJ then a is arithmetic, meaning that cr is a 
discrete measure supported by a geometric sequence. Note that in the conservative case (3* = 1 and in 
the dissipative case (3* <1. 

Equation Q has no real solutions (thus the Malthusian exponent is not defined) only if (/)(/3a+) < 1. 
An example of this situation is the measure 

a(Ax) — c i{x<i/2} x'^^l'^ log~^ X Ax 
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with a suitable choice of c, and (5a = 1/2. 

Further assumptions about the reproduction law will be introduced in a due place. Specifically, the 
i^-convergence results in Sections 12.21 and 14.11 require that 



Example. Consider a dissipative reproduction law induced by the uniform stick-breaking. Let 
C/qj C^ij • ■ • be i.i.d. uniform, and let {CjjjeN stand for the rearrangement in the decreasing order of the 
sequence 

{l-U,)\[Uk, J = 1,2,... 

meaning that the uniform portion of the size 1 — C/q is lost, and the rest is fragmented by the 'random 
alms' principle (as Halmos called stick-breaking). Equivalently, the offspring {'fj}jeN of a unit particle 
can be seen as the collection of sizes of a random Poisson-Dirichlet distribution with parameter 1, upon 
removing a size selected by a size-biased pick. Building the characteristic function 



oo _^ 



^^(l + /?p-+i /?(/3 + l) 
we see that the abscissa is at /3a = and the Malthusian exponent is (3* — (—1 + a/5)/2. 



It should be noted that there is no substantial constraint on a imposed by the requirement that a 
be a structural measure. Given a on [0, 1], satisfying cr]0, 1] > 1 and <t{1} < 1, a possible reproduction 
law satisfying ^ with this structural measure can be constructed as follows. Decompose cr = cti + (T2 so 
that (Ti be probability measure and a2 some other measure. Let be a random point with distribution 
CTi and let ^2; • ■ ■ be the atoms of a Poisson point process on the unit interval, with intensity measure 
(T2. Clearly the point process will have intensity a. 



2.2 Genealogical structure and the intrinsic martingale 

In this section, we develop some elements on the genealogical structure of the fragmentation. This can 
be viewed as a different parametrisation of the process, in which the natural time is replaced by the 
generation of the different particles. Specifically, we consider the infinite tree 

oo 

U -.^ U N" , 

n=0 

with the convention N'^ = {0}. The elements of U are called nodes. For each u = (ui, . . . , u„) G we 
call n the generation of u and write |u| = n, with the obvious convention |0| = 0. When n > 1 and 
u — (ui, . . . , u„), we call u— := (wi, . . . , Un-i) the father of u, and u,i = . . . , u„, i) the i-th child of 
u. 

It will be assumed that a > and that the fragmentation process X starts from a single particle with 
unit size, unless explicitly indicated. We encode X by putting marks on the nodes of the infinite tree K 
as follows. The initial particle with unit size corresponds to the ancestor 0, and the mark of is the triple 
{i$,g$,d(ij) = (1,0, e) where e is the instant of the first reproduction, so that [50,d0[ is the time-interval 
during which the ancestor particle is alive. The nodes of the tree at the first generation are used as the 
labels of the particles arising at the first split, i.e. ^i, . . . . . .. Again, the mark associated to each of 

these nodes z G N^, is the triple (^i, gi,di), where gi = stands for the birth-time of the i-th child of the 
ancestor, and di = gi + ^fe^ for its death-time. And we iterate the same construction with each particle 
at each generation. 

Plainly, the description of the dynamics of fragmentation entails that its representation as a random 
marked tree enjoys the branching property. Specifically, the distribution of the random marks can be 
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described recursively as follows: Given the marks {{£,v , gv , dv) , \ v\ < n) at nodes of the first n generations, 
the marks at nodes of generation n + 1 can be expressed in the form 

iiu,gu,du) = iiu-iu, du~ : du- + ^u^u) , u = (ui, . . . ,ti„+l), 

where u— = (ui, . . . , u„) is the father of the node u, and 

• when u— = {ui, . . . , u„) G N" describes the nodes at the n-th generation, ($ui,...,tt„,i, « G N) are i.i.d. 
random sequences distributed according to the reproduction law, 

• when u describes the nodes at the (n+ l)-th generation, the variables are i.i.d. standard exponential 
variables which are independent of the sequences - 

One says that extinction occurs when £,u — for all nodes u at generation n for some large enough 
n E N. The assumptions ensure that the process Z„ := #{u G N" : ^„ > 0} is a supercritical 
Galton- Watson process (possibly taking the value oo), so the probability of non-extinction is strictly 
positive. 

It is well-known from the works of Jagers |23, Nerman (31j and many other authors, that the Malthu- 
sian hypothesis for branching processes is connected to a remarkable martingale. The latter is often 
referred to as the intrinsic martingale, it plays a crucial role in the analysis of the asymptotic behaviour 
of branching processes. The following statement, expressed in terms of the genealogical coding, is part 
of the folklore (we refer to for the last part of the claim). Recall that (3* denotes the Malthusian 
exponent. 

Proposition 1 Under assumptions ^ and |j2Jl, the process 

\u\—n 

is a martingale which is bounded in and in particular, uniformly integrable. Its terminal value Moo 'is 
strictly positive conditionally on non- extinction, and satisfies the distributional identity 

oo 

Moo = Y.€^^^ (4) 

i=i 

where are independent copies of Moo, also independent o/{^j}jgN- The identity taken together with 

conditions (O and EA/qo = 1 characterises Moo uniquely. 

Remark. Observe that in the important case when the reproduction law is conservative, in the sense 
that — ^ 21. s., we have /3* — 1 and A/„ — 1 for all n G N, so that the statement is trivial. 

So far we have described the fragmentation process in terms of its genealogy; however the problems 
of interest are often expressed in terms of the natural time scale. More precisely, the configuration of the 
fragmentation process at time t > consists in the set {Xj{t)}j^fq of particles coexisting at time t > 0, 
that is in terms of random point measures 

Recall we assume that the process starts with a sole particle of unit size, that is X{0) = {1,0,...}. 
Denoting X^^^ the fragmentation process that starts with a particle of size j/ > 0, we have the fundamental 
self-similarity identity 

X^y\t)^yX{ty'^). (5) 

The intrinsic martingale is indexed by the generations of the infinite tree. In terms of the time scale 
of the fragmentation process, we write {Tt)t>o for the natural filtration generated by the fragmentation 
process X = {X{t),t > 0). Proposition ^ then yields the following result which could also be proved by 
techniques of so-called stopping lines applied to the branching process induced by the genealogical coding 
(see e.g. Theorem 6.3 and Corollary 6.6 in Jagers |21J). 
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Corollary 2 Under the assumptions of Proposition^ for every t>0, we have 

E (Afoo I ^t) = M{t, /?*) J2 ^/'(*) ' 

J 

hence the process M{-,f3*) is a square- integrahle {J^t)-fnartingale with terminal value Moo- 
Proof. We know that M„ converges in L^(P) to Moo as n tends to oo, so 

^{Moo I ^t) = lim E(M„ I Tt) . 
On the other hand, it is easy to deduce from the Markov property apphed at time t that 

oo 

E(M„ I = (t)i{G(x,(t))<„} + ^ri{<i.<t} ' 

i=l |ti|=n 

where G(x) stands for the generation of the particle x (i.e. G(^u) = \u\), and d„ for the instant when the 
particle corresponding to the vertex u splits. However, for each fixed vertex u d U, du is bounded from 
below by the sum of |u| + 1 independent exponential variables which are independent of ^u- It follows 
that 

lim E J2 l{rf„<t} =0> 

\u\—n 

and we conclude that E {Moo \ ^t) = M (t, /?*). □ 



3 Asymptotics in mean 
3.1 Power sums and their means 

In the sequel, we shall find useful to consider the power-sum functionals 

M(t,/3)^^xf(0 
j 

and their means 

m{t,l3) ^EM{t,l3). 

For shorthand we sometimes refer to the size of a particle raised to the power (3 as the l3-size, thus 
M{t, (3) is the total /3-size of the population existing at time t. Two instances with obvious physical 
interpretations are the 0-size equal to the number of particles, and the 1-size equal to the total size of 
the ensemble. 

As the set of particles alive at time t is a part of the set of particles that are born before time there 
is the inequality 

oo 

m{t,P) = E^ Xf (t) < E ^ ef l{s„<t} , 

i=l u&A 

where gu denotes the birth-time of the particle labelled by the node u. Observe that for each node u 
at generation n, g{u) can be bounded from below by the sum of n independent standard exponential 
variables which are also independent of and that, in the notation of Section l2.ll 

\u\—n 

It now follows from classical large deviation estimates that m(t, P) < oo for every t > whenever (3 > j3a- 
The first-split decomposition with application of shows that M(t, (3) satisfies the distributional 
identity 

M{t,l3) ^ + l{,>rf,} M,{^'^{t - dg)) (6) 

3 
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where d0 is the exponential life-time of the progenitor, and the Mj's are independent replicas of M(-, f3), 
which are also independent of and {^j}. Computing expectations we arrive at the integral equation 



Jo Jo 

Differentiating we see that m(-, /?) is a solution to the Cauchy problem for the integro-differcntial equation 

dtm{t,(3) = -m{t,f5)+ [ mix^t, f3)x'^ a{dx) . (7) 



which must be complemented by the initial value to(0, /3) = 1 . Uniqueness of C°° solutions for equations 
of this type is shown in 20 . 

Equation ||7J| defines functions m for all > Pa- For the higher derivatives we have 

d^m{t,f3) =d^m{0,f3) m{t,kp + a), 

thus m is increasing in t for (3 < (3* and decreasing for (3 > (3* . 
Solving |7J in power series is straightforward. Introducing 

i^{i3) = i- m 

and then 

n-l 

l{n,(3)=l[iji(3 + ak) (8) 

k=0 

(by convention, 7(0, /3) = 1), we compute 

m{t,fi)^Y.-^l{n.P) (9) 

n=0 

which is an entire function of t G C. It is indeed the right solution because from the formula for derivatives 
it is clear that m(-, (3) should be C°° for t > 0. See for yet another interesting representation of m(t, f3), 
in the form of a generalised Dirichlet series. 

3.2 A contour integral 

When (j){/3) is a rational function, splitting 7(- , (3) in linear factors shows that the series © represents 
a generalised hypergeometric function, in which case the \t\ — > 00 asymptotic expansions have been 
thoroughly studied by Mellin-Barnes' contour integral technique, see [23 ■ This method extends to the 
more general situation considered here (also see 

Call P singular if ?/;(/3 + <^n) = for some integer n > 0. For singular /? the series m{t,j3) is 
a polynomial, thus 171(1, P) = j{n,f3)t^ + 0{t^~^) for t 00. Analogous asymptotics hold also for 
nonsingular (3, but it is more difficult to justify, because m is then an infinite series which starts alternating 
from some term. A good heuristic amounts to substituting m ~ ct°' into Q) - the left-hand side is then 
of the order t°-^^ while the right-hand side is o(t°) exactly when ip{/3 + aa) — 0, which suggests that 
a — {[3* — (3) /a is the right exponent. Although this kind of reasoning can be made precise it gives no 
idea of the coefficient, see |2iJJ. 

Assuming (3 nonsingular we extrapolate the function 7(- , [3) from the integer values to arbitrary 
complex values z (such that (j}{az -f (3) is defined) by means of the formula 

Convergence of the product follows as in Section 5. Thus defined, 7 satisfies the functional equation 

7(z + l,/3) = V(/3 + a2)7(z,/3) 
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reminiscent of the well-known equation for Euler's gamma function. 

All singularities of the function 7(- , /3) are the poles located at roots of ||2Jl. Let 



Vp = {z: 3n>0, tP{/3 + a{n + z)) = 0} 

be the set of singular points. Because has no solutions to the right of /?*, the rightmost point of Vp 
is z/s := {13* — /3)/a, where 7(- , fi) has a simple pole provided (3* > Pa- 

Still assuming /3 nonsingular we have Vp n {0, 1, . . .} = 0. Since the poles of r(— z) are nonneg- 
ative integers, the function r(— 2)7(2, (3)V' (with t as parameter) also has these poles, with residue 
(— l)"<:"'7(n, /3)/n! aX, z = n. Defining C to be a vertical line between and np := min(0, [5Rz^]) 
we obtain by the residue theorem and an estimate of 7 

tlL ^{n, P) = ^ I r(-z)7(z, I3)z' dz . (11) 
^-^ n\ 27ri Jr 

n=ni3 

If (3* > (3a the function ,/3) is meromorphic in an open strip containing the line 5Rz = zp, and the 
residue at zp is 

as it follows from the identity 

ip{(3 + as) 

upon expanding the ratio. Replacing C by another integration contour C located in the half-plane 
3?z < ^zp so that all poles of 7(- , j3) in this half-plane lie to the left of C we obtain the principal-term 
asymptotics of to. 

Theorem 3 Suppose (3* > (3a, the structural measure <j is nonarithmetic and conditions hold, then 
\ a J ayy[(j*) \ a J 

for ^(3>pa- 



Res^^7(z, P) = 7 ( — I — , a + P 



Using the identity j{—z , az + P) 7(2, /?) = 1 we can re- write the 7-factor in H12|) as 

1 

j{{P-P*)/a,a + P*) ■ 

The restriction of "0 to the real segment ]Pa, oo[ is plainly a concave increasing function, so the 
condition P* > Pa entails < i{j'{P*) < 00. We also remark that in the arithmetic case other poles 
on the line 5Rz = /?* would contribute to the coefficient. Further terms of the asymptotic expansion 
can be obtained by pulling the integration contour through other poles left of P* , as long as admits 
a meromorphic continuation, which can go beyond the convergence abscissa. If P* is the only point of 
Vp in a closed strip P* ~ < ^z < P* then the rest-term in (O is estimated as 0(t(/3*-/3)/a-e) ^j^j^ 
e — min(l, 6/a). 

3.3 Convergence of the mean measures 

We encode the configuration of sizes X{t) = {Xj{t)} into the random measure 

3 
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The associated mean measure cr* is defined by the formula 

/ f{x)a:{dx)=Ej2f{t'/"X,{t))xf{t) (13) 
Jo ^ 

which is required to hold for all compactly supported continuous functions /. It is easily seen that is 
a probability measure. Our next goal is to show that the measures CTj converge weakly to a probability 
measure p on ]0, cx)[ . 
Because 

POO 

the convergence of t'''"^ '>^°'m{t, (3) implied by (|12|l amounts to the convergence of power moments 

Jo ' V « y a^'(/3*) 7 ((/3 -/?*)/", " + 

Specialising f3 — (3* + ak this simplifies to 



in application of It is easy to check and is well known that for n = 0, 1, . . . the related 

moment problem is determinate, whence the following result. 

Theorem 4 Under assumptions of Theorem |31 the measures CTj converge weakly, as t ^ oo, to a 
probability measure p uniquely determined by its power moments 



Theorems 13 and 01 imply that 

til^-t^n/o^jnit^P) ^ I xl^-^' p{dx) (15) 



which extends the convergence of expectations in (|13|1 to a wider class of functions /. 
3.4 Self-similar stick-breaking process 

A key tool in the conservative case treated in ^| I18L [SI E| |S| has been the following observation related, 
somewhat paradoxically, to the simplest dissipative case of a singleton ensemble with reproduction law 
{rj}, where 77 is a random variable assuming values in ]0, 1]. That is to say, if at some time the particle 
has size x then, independently of the history, the particle shrinks with probability rate x" and the new 
size after the shrink becomes xrj where rj follows a, and a is the probability law of 77. 
We recollect briefly a result from 18 (also see [51 E]). Introduce 

^(/3) = 1 - / xf^ a{dx) 
Jo 

and suppose "0(0+) < 00. Let Lt be the sole size at time t, and fh{t,j3) — Eif. Assuming a non- 
arithmetic, Brennan and Durrett |18| proved that for t ^ 00 

t^/^Lt ^ Y^'" and m(t, /3) -> E F'^/" 
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where y is a random variable with moments 



The convergence of moments m(t, (3) was shown for all real [3 strictly to the right of the convergence 
abscissa of (which is nonpositive due to the normalisation (t]0, 1] = 1). They also suggested the 
explicit representation 



oo 



k 



^=E^^-n^" (16) 

fc=0 j=0 

where all 7]j, ej are independent, ej is mean one exponential, rjj for j > are replicas of rj, and rjo follows 

the law ^ 

™/ , N o^\x, lldx 

Each product in (|16|l corresponds to the size of the particle after k splits, conditionally given the initial 
size is 770, thus formula H16|) identifies Y with the well-known exponential functional of a (stationary) 
compound Poisson process, see e.g. the survey |13| . 

In the conservative fragmentation case, the above 'stick-breaking' process describes the evolution of a 
particle tagged by an atom of isotope that was injected at a random uniform location into the progenitor 
unit size. The mechanism which determines the line of descent of the tagged particle amounts, at each 
consecutive split, to a random size-biased pick from the child particles. Thus, defining a{dx) :— x(T{dx) 
to be the distribution of a size-biased pick from {Cj}, we obtain the relation m{t, (3 — 1) ~ m{t, j3) which 
was observed in J_8,, p. 112. 

In the general nonconservative case choosing a{Ax) :— x^ a{dx), we still get fn{t,(3 — (3*) = m{t,j3) 
and = "ipiz + An interpretation of these relations akin to the tagged fragment process is 

available through the so-called spine which appears in the conceptual approach of Lyons et al. |26II25| for 
convergence of martingales in branching processes. The measure p may be identified with the distribution 
consequence of this discussion and the result of Brennan and Durrett is the following corollary. 

Corollary 5 The conclusion 0/ Theorem ^ remams valid even if the assumption /?* > Pa is replaced 
by the weaker ^p'{f3*+) < 00 



The tiny improvement upon Theorem^] appears in the case where the characteristic function is defined 
in a closed half- plane and [3* — (3a, i-e. the Malthusian exponent falls exactly on the convergence abscissa. 
An example of such situation is the structural measure of the form a{dx) = cl^x<i/2}X~^^'^\ logxl""^ da; 
with a suitable c. The method based on contour integration requires in such cases the analytical contin- 
uation of in a domain to the left of f3a- 

Alternatively, along the lines in jl7l 1181 IIT)] , the renewal theory can be applied also in the noncon- 
servative case, to prove first the convergence of measures CTj — > p and then to justify the asymptotics of 
mean power-sums using the uniform integrability. 

4 Limit theorem for the empirical measure of the fragments 
4.1 L^-convergence 

Our principal result improves on the convergence of the mean measures Ct — » p in Theorem 0] and says 
that the scaled empirical measures induced by X(t) converge in a L^-sense to the measure M^oP, where 
Moo is the terminal value of the intrinsic martingale (cf. Proposition 
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Theorem 6 Assume 0, that /?* > (3a and that a is nonarithmetic. Then for any bounded 
continuous f 

L^-lim V Xf {t)f{t^'''X,{t)) = Moo / f{x)p{dx). 
Proof. We need to show that 

^ ^ '"^ f{x)p{dx)\( r g{x)p{Ax)] (17) 



for positive / and g bounded from above by 1. Indeed, suppose l|17|l is shown. Denote 



Take f — g to conclude from (|17l) that 

EA? = EAfi ( I f(x)p(dx) ^ 



hm EA^ = EAf^ ( / fix)p{dx) 
Similarly, setting g — 1 



lim E{AtM{t, 13*)) = EAf4 / fix)p{dx) 



t — >oo 







Recalling from Corollary |2 that EAf(t, /3*)^ — > ^M"^ and combining the above we get the desired 

lim E (At - M{t, (3*) I f{x)p(dx)\ = . 



To prove (|17|l let us replace t hy t + s and condition on the configuration of sizes X{s). At time 
i + s two coexisting particles may stem from the same ancestor that lived at time s or from two different 
ancestors; write i j in the first case, and write i j in the second. The sum in the left-hand side of 
(|17|l is split then in two 

S-i + ^2 = E I E • • • \X{s) I + E I 51 ■ ■ • 

Using the fundamental self-similarity relation 10) and the Markov nature of the fragmentation process 
we estimate the first sum as 

k \ 1 

hence by (fT^ and Corollary 12 

ESi < const s"''*/" ^0 as s ^ 00 

uniformly in t. 

Dealing with S2 requires more effort. We use the parallel notation yj — Xj{s). Write i \ k if the 
size Xi{t + s) stems from yk- By independence, the descendants of different particles with sizes yt and 
yi evolve independently, thus grouping the sizes Xj(t + s) by the ancestors at time s yields 
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However, by self-similarity and convergence of the mean measures 

.V I, "'0 



^Y.yT iiyDaiii + sY'^y,X,(ty1)) ^ yf / g{x)p{Ax) 
as i — > oo, therefore by dominated convergence 

ES2^ (^J^^ f{x)p{dx)^ (^j^ g{x)p{dx)^^Y.xf[ {s)xf {s) 



as s — > oo. It remains to note that 



E ^ * {s)xf (s) ^ E 5] X[ {s)xf {s) ^ em! - 

kjte k , I 

because 

E5]xf*(s) = m(i,2/3*)^0. 



□ 



Remarks. We mention that if we replace the assumption Q by the weaker 



for some 1 < p < 2, the calculation of the expectation of the sum of the p-th powers of jumps then 
shows that the intrinsic martingale M„ in Proposition ^ is bounded in L^, see e.g. Neveu [321. Then 
techniques of Nerman based on generalisations of the law of large numbers for branching processes, 
can be adapted to extend Theorem to this situation (of course, Li^ convergence then has to be replaced 
by L^'-convergence) . 

Finally, it is known that in the binary conservative case the scaled empirical measures converge with 
probability one ^Hl- It would be interesting to extend this result to the nonconservative case. 

In parallel to p5|l there is the following extension of Theorem |S1 to power functions. 

Corollary 7 Under assumptions o/Theorem|Hl 

L2- hm t^^-f^'^'^y X^{t) = Moo / xP-^'p{dx) 

Proof. Along the same line, the proof is reduced to showing that 

is bounded uniformly in y G ]0, 1[ . And the latter follows by noting that ((SJ and (|15(1 imply 

, . Q poo 

t(/3-r)/«EV(xj^)(i))^/* / x^-f^'pidx). 

□ 
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4.2 Comparison with Branching Random Walks in continuous time 

In this section, we compare the case of positive indices of self-similarity a > which we have considered 
so far, with the simpler homogeneous case when a = 0. 

So suppose here that a = 0. The process Z{t) — {Zj{t)} with Zj{t) ~ — log Xj{t) is then a continuous- 
time branching random walk, as studied in |35II16| . We have in the notation of Section [3. II that "f{n, j3) — 
V'"(/3) and the identity © for the moments of power sums thus becomes m{t,(3) — exp(— a 
formula which is well-known in the context of branching random walk; see 

It is read from the work of Biggins |15[I16| that for every p > (3a, the process 

Wit, P) := exp(t^(/3)) ^ , t>0 



is a martingale with cadlag paths, which converges almost surely and in their mean as t — > cx3 provided 
that i^iP) < Ptp'iP), see JT] for details. Note that this holds in the special case when f3 — P* is the 
Malthusian exponent, for which there is the identity W{t,P*) ~ M{t,P*) ~ Mt. 

When furthermore the reproduction law is not arithmetic, the asymptotic behaviour of the empirical 
distribution of sizes can be described as follows: for every continuous function / : R — > M with compact 
support. 



hm Vte-*('5^'(^)-'^(^)) V f{ti:'(p) - Z,{t)) - 1^^°°'^^ r f{-z)e^^ 



j 



Mz. (18) 



where W{oo,P) is the terminal value of W{t,P). See 

Theorem |B1 bears obvious similarities with H18() . It is interesting to observe that in the homogeneous 
case a — 0, sizes decay exponentially fast and the limiting scaled empirical measure is always exponential 
(up-to a random factor), whereas for a > the decay of sizes is power- like and the limiting scaled 
empirical measure depends crucially on the structural measure a (more precisely, a can be recovered 
from the limiting scaled empirical measure for a > 0, but not for a = 0). 

5 Examples 

5.1 Fihppov's example revisited 

Extending an example in Filippov (see |1H| . section 8) consider the structural measure 

a{dx) ^ Xx"-^ dx , j:e]0, 1] 

with parameters A > min(0, 0) and arbitrary 6 Cz R. For X < 9 + 1 the fragmentation is mean-value 
dissipative. We have 

The characteristic function is thus meromorphic in C with a unique simple pole at Pa — —9. 
Computing 

we see that this is the ratio of two Pochhammer factorials, with A = {P ~ P*)/a and B = {9 + P)/a, 
thus m(t, P) ~ iFi{A-, B; —t) is Kummcr's hypergeometric function. The analytical extension of 7 is 

Computing the moments we obtain 

00 

p{dx) = (A/a)fc 
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which identifies p as 



pidx) = . dx, a;>0. 

1 (A/a) 



Note that the shape parameter 9 cancels and does not appear in p. It foUows that at := x ^ al (the 
intensity of M Sii/aXj{t)) converges to the measure 

x-^^p(dx) = — a;^-ie-^° dx, x>0 
r(A/a) 

in accord with the case A = 2, = 1 considered in |18j in connection with the conservative binary 
fragmentation with uniform and ^2 = 1 — ■ 

It foUows that the mean number of particles satisfies 

r(A/a) 

which agrees with a special case in ^^l- Of course, this formula makes sense only for 9 > 0, because our 
asymptotics for m{t, /3) hold only for 5ft/3 > /3a, thus for 9 < the value P = is not considered. 



The case A = 1, 6* = 0, when a{dx) — x~^dx corresponds to the conservative fragmentation generated 
by the uniform stick-breaking, as in the example in Section [2. II (but without removing a piece). It is well 
known that the distribution of a size-biased pick from {^j } is uniform, and this implies that the intensity 
measure of ^ (5^^ is indeed a(dx) = x~^dx. 

5.2 Hypergeometrics 

The following is a further generalisation of Filippov's example, and covers the class of dissipative binary 
fragmentations treated in 0. Consider a Dirichlet polynomial 

5(a:)=^A,x^^-i (19) 
which is non- negative on ]0, 1] and has real parameters satisfying 

j=i 3 

Then a{dx) — g{x)dx is a measure on ]0, 1] with rational characteristic function 

and by the assumption the rightmost root of (f>{P) = 1 is positive, denote it also /3i = /3* and denote 
further roots (32, ■ ■ ■ , Pp (the roots are certainly different from the poles of <j>). 
Observe that 



1 /3 + 

thus assuming a — 1 (without loss of generality) we have 



n=0 j- 



J] jn 
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where we recognise a generalised hypergeometric function of the type pFp. By Q we have mlt, 13) 
c{l3)t'^-P' for > /3a. Noting that 



/3* - 0j 

13* + 01 H /3* + e, 



J=2 

and manipulating infinite products, the coefficient is evaluated in terms of the gamma function as 
This allows to recover the density by Mellin inversion as 



x>0 



which is an instance of Meijer's G- function, see |27| . 

The limit measure is uniquely determined by the integer moments which can be written as 

(fc-1)! /r {(3* + i + 0j)k-i 



where the derivative may be also computed as 



6 Fragmentations with infinite reproduction measure 

We sketch how the preceding results can be extended to a class of self-similar conservative or dissipative 
fragmentations with infinite reproduction measure. Such processes were introduced in j2) where the 
reproduction law was called 'dislocation measure'. 
Let V he & measure on the infinite simplex 

C30 

A = {(s,) : > 0, iO, ^s, <1} 

such that the integral 

m--^ I (i-E^.^)Kds) (20) 

satisfies 1 < il){(3) < oo for some /3 > 0. For a > we can define a fragmentation process {X{t),t > 0) 
with the property that a generic particle of size x gives birth to a collection of particles of sizes xsj with 
s ~ {sj) G B at rate x°^v{B), where Sj runs over nonzero coordinates of s and i? C A runs over Borel 
sets of finite i/-measure. 

If is a probability measure it can be regarded as a reproduction law by defining (^j ) to be a random 
element of A with distribution v. In this case the structural measure is identified with the superposition of 
marginal distributions of and the definition (|2(J|I agrees with our definition of the characteristic function 
in Section im The case < oo is easily reduced to the case i^(A) = 1 by the obvious time-change. 

In the case ^'(A) = oo some features of the fragmentation process are different, in particular, each 
particle produces infinitely many generations within arbitrarily small time period. As a consequence, 
the life-time of a particle is not a well-defined quantity, we do not have a tree representation for the 
genealogical structure, nor equation like Still, we can define j3a and (3* exactly as in the case of finite 
measure, and consider /3-sizes for 3? /3 > /Jo- Formula Q remains valid and can be proved by an argument 



14 



exploiting approximation of v by suitable finite measures, or by using the methods developed in jl()lll2| . 
For a > conclusions of Theorems 0| and |^ remain valid if we assume that 



< cx) and i \^sf -l \ v{ds) < oo , 

and impose a non-arithmeticity condition on v. The analog of Corollary |2] is shown by arguments similar 
to those in Theorem 2 of |S] , and the analog of distributional equation generalises in the form of the 
Laplace transform identity 

I \{C{esf )v{ds) (21) 

for C{9) := Ecxp(— 0Moo)- The question about the uniqueness of solution to (|21|l with infinite measure 
seems to have not been considered before and remains open. 

Finally, we observe the obvious interchangeability between parameters a and f3 (which holds no 
matter whether the reproduction measure is finite or infinite) . In the variables = the fragmentation 
process has the life-time parameter a/ (3 and differs only by a particle-wise transformation of sizes. This 
observation enables us to handle situations when the measure v is not supported by the infinite simplex 
A, but rather by 

oo 

A^ = {(s,):s, >0, s,iO,J2^j<^} 

j=i 

for some f3 > 0. Note also that the change of variables with /3 < yields mathematically equivalent 
(though physically curious) process where the individual 'fragments' grow, but the decaying life-times 
slow down the total increase of size. 
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